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The scaling behavior of cyclical surface growth (e.g. deposition/desorption) , with the number 
of cycles n, is investigated. The roughness of surfaces grown by two linear primary processes 
follows a scaling behavior with asymptotic exponents inherited from the dominant process while the 
effective amplitudes are determined by both. Relevant non-linear effects in the primary processes 
may remain so or be rendered irrelevant. Numerical simulations for several pairs of generic primary 
processes confirm these conclusions. Experimental results for the surface roughness during cyclical 
electrodeposition/dissolution of silver show a power-law dependence on n, consistent with the scaling 
description. 
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Much interest has been devoted in recent years to scaling phenomena in kinetic growth of self-affine surfaces. They 
are observed in a variety of phenomena such as crystal growth, vapor deposition, molecular-beam epitaxy, electro- 
chemical processes, bacterial growth, burning fronts, etc. H p] 0]. Similar rough surfaces may also be generated 
^ by reciprocal processes of surface withdrawal caused by desorption, corrosion, evaporation, dissolution, and alike [pi. 
In many natural and artificial systems of interest, however, surfaces are not formed by a single process of growth or a 
sole process of recession. Rather, they are the product of a combination of both. We investigate here the fundamental 
scaling properties of surfaces formed by cyclical growth processes. We focus on surfaces formed by two alternating 
primary processes. Growth/recession cycles (excluding the trivial cases in which the primary processes are time- 
reversed images of each other) are our main interest but cycles of two different growth (growth/growth) processes will 
be addressed as well. Cyclical behavior, prevailing in all natural phenomena, may be found in many of the systems 
cited above. Typical examples include weather and light affected processes in organic (as the expansion/curtailment 
of a grass lawn or a bacterial colony according to the availability of water or nutrient) or non-organic (like alternating 
underwater erosion/sedimentation whether or not the water is flowing) systems. They are also widespread in techno- 
-y-s logical applications (such as rechargeable batteries for which shortcircuit by the metal accumulated on the electrodes 

,_^ ' is one of the failure mechanisms). Understanding the cyclical scaling properties may lead to accelerated testing, and 
t-H , performance improvement, of such systems. 

Our main challenge is to generalize the scaling approach in order to make it applicable to cyclical growth processes. 
The analytical, numerical, and experimental investigations, summarized below, lead us to the conclusion that this 
is indeed possible, provided that the number of cycles n substitutes for the time variable t in the scaling relations. 
We explore this generalized scaling behavior in several cyclical processes and reach some general conclusions on their 
behavior and their relations to the scaling relations of the primary processes, of which they are composed. 

In the primary processes, the surface width W(L,t), where L is the lateral size of the system, is defined as 
W(L,t) —< (h(r,t) — h(t)) 2 >2. In this definition h(r,t) is the surface height and h(t) —< h(t) >— vt is the average 
height, with v being the average growth velocity. W(L,t) obeys the following scaling form M M: 
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W(L,t)~L a g(L/at)), (1) 



t;(t) ~ i 1 / 2 is the lateral correlation length. For large time t » L z : W ~ L a , while for t « L z : W ~ t 13 , where 
/3 = a/z is the growth exponent. A related relation holds for the mean square height difference < (h(r + x,t) — 
j_j ■ h(x,t)) 2 >= 2(C(0,t) — C(r,t)), where C(r,t) is the (equal-time) height-height correlation function. 

The growth processes fall into different universality classes |2) . All processes within one class share the same expo- 
nents and their asymptotic continuum stochastic equations differ at most by irrelevant terms (in the renormalization 
group (RG) sense). Using the symbolic index i = 1, 2, .. to denote different processes, the ones we consider here follow 
growth equations of the form: 

dH dt^ =Mh}+Vi(r,t)+Vi, (2) 

where Ai{h} is a local functional depending on the spatial derivatives of h(r,t) and < rji(r,t)r]i(f / ,t') >= 2DiS(f — 
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We recall the simplest generic growth processes of this type: (i) Random deposition (RD) for which Ard = and 
f3 = 1/2 (a and z are not defined); (ii) The Edwards- Wilkinson (EW) model J7| of preferred growth at local minima 
has A E w = ^V 2 /i and a = ^=^, (3 = ^p, z = 2; (hi) The Kardar-Parisi-Zhang (KPZ) M model which accounts for 
the growth being locally normal to the surface has Akpz = v\J 2 h + 4(v/i) 2 and a — 1/2, (3 — 1/3 and z — 3/2 in 
2d, while a ~ 0.39, (3 ~ 0.24 and z ~ 1.61 in 3d; (iv) The Das Sarma-Tamborenea (DT) g (see also fnj) for MBE 
deposition or growth on kink sites has Aot — v^J^h with a = ^2, /3 = ^p, and 2 = 4. 

We investigate theoretically and experimentally the hypothesis that in cyclical processes the scaling law in Eq. (1) 
should be replaced with: 

W c (L,n) ~ L a g c (L/Un)) (3) 

with £ c (n) ~ n 1 / 2 . 

We begin with the study of cyclical processes composed of two primary linear processes, namely those for which 
Ai{h} — ai(\)h(f, t), where Oj(V) is a differential operator (e.g. the EW and DT models) and for which time-reversal 
symmetry holds if the height is defined w.r. to the average height (i.e. for h(r,t) — h(t)). The first process in the 
cycle (i = 1) is of duration T\ = pT and second process (i = 2) lasts T 2 = (1 — p)T. The total duration of one cycle 
isT = T 1 +T 2 . 

We also define f(t) as the fractional part of t/T. The growth equation thus becomes: 

dh 

— = [a x h + mMp - /(*)) + [a 2 h + m ]e(f(t) - p), (4) 

where Q(x) is the unit step function. 

For such linear processes, the full scaling behavior may be retrieved by looking at h(q,t), the Fourier transform 
(FT) of h(r,t). The first observation we make is that in these time-reversible processes, only the averaged height 
h(t) = h(q = 0, t) is sensitive to the difference between growth (vi > 0) and recession (vi < 0). In terms of the average 
velocity v c = pv\ + (1 — p)v2, it is given by: 

h(t)/T = nv c + v 1 f(t)Q(p - f(t)) + [(vi - v 2 )p + v 2 (f(t)]e(f(t) - p). 

The roughness is insensitive to the sign of Vi and hence will not discern between a growth/growth and a 
growth/recession cyclical process (as long as a, and Di are not altered) . 

In Fourier space, the growth equations for the modes with q ^ may be integrated. The h(q,t) after n (res. 
(n + p)) cycles are assigned as the initial conditions for the (n + l) th application of the first (res. second) process. 
The structure factor S(q, t) =< h(q, t)h(—q, t) > (FT of C(r, t), for simplicity we assume spatial isotropy in the basal 
plane) is then derived by averaging over the noise. We define di = Oj(g)Tj and a c = [a±p + 02(1 — p)]T, in terms of 
which we find S c (q,n) = S(q,nT) after exactly n cycles to be: 

S c (q,n) =exp{-2a c n}S(g,0)+ [^ exp(-2o 2 )(l - exp(-2oi)) + ^(1 - exp(-2o 2 ))] [ ^jg^g ] . (5) 

A similar expression may be obtained for S c (q,n + p) and straightforwardly extended to any time t = (n + f)T. 
For small q, s.t. a c (q) << 1, S c (q,n) takes the form: 

S c (q, n) - -^ [1 - exp(-2( ac (<?)Tn)] , (6) 

a c (q) 

where we introduce the effective parameters of the cyclical growth process: D c = pD\ + (1 — p)D 2 , and a c (q) = 
pa\(q) + (1 — p)a 2 (q) = a c /T. The same coclusion may be reached from coarse graining Eq.(4) and eliminating all 
modes with frequencies larger than 2n/T |11[ (this also yields the corresponding propagator G c (q, n) on time scales 
larger than T). In terms of these effective parameters, S c (q, n) (G c (q, n)), with large n = t/T, of the cyclical process 
are equivalent to S(q,t) (G(q,t)) of a primary linear process. Hence, the scaling behavior, presumed for the former 
in Eq.(3), indeed substitutes for that of Eq.(l) which holds for the latter M. 

In the n — > 00 limit, the asymptotic large-sacle roughness is determined by the small q divergence of S c (q, t). Since 
a i(<l) ~ \l\ Zi i h is the process with the smaller z% which dominates the large L cyclical roughness: W C (L, 00) = AL a , 
with a — min(ai,a 2 ) = min(zi,z 2 ) — (d — 1) (the amplitude A is proportional to D c and is determined by both 
primary processes). The larger on appears as a correction to scaling exponent (whether or not it is the leading one 
depends on how its contribution compares with that of a potential subleading term in the dominating dj ((?)). Note 
that a subleading term in a c (q) might affect the behavior on a smaller scale, if its amplitude is large. In that case, the 
leading behavior takes over only beyond a crossover length (at which both contributions are comparable). As could 
be expected, the longer the non-dominant process lasts, the larger is the crossover scale to the dominant behavior. 



As for the dynamic exponent z of the cyclical process, since n appears always multiplied by a c (q), the slower of 
the two primary processes dictates the cyclical dynamics: z = min(zi, 22)- The correlation length increases with n 
as: £ c = n l / z with an amplitude proportional to [X^] 1 / 2 , with Ti of the dominating process. For the initial cycles 
(nT << L z ), the roughness grows as nP , with (3 = a/z = (z — (d — l))/z. 

For processes described by stochastic non-linear equations, the asymptotic behavior is explored by the RG approach 
||. An approximate RG procedure for cyclical processes may consist in first coarse graining the free cyclical 
propagator (obtained from the two primary free propagators as explained above) until it becomes that of an effective 
linear, non-cyclical process. Using this effective linear process as the "free" part, all the non-linear terms of the 
primary processes, with bare couplings multiplied by p (or 1 — p), are added as "interactions". The RG flows may 
then be derived following the ususal steps pT| . In this approximation the initial flow of the couplings is replaced by a 
simplified one. It is implicitly assumed that simplfying the initial flow will not alter the ultimate fixed point for each 
of the renormalized couplings. Although this seems very plausible, it might not hold for all cyclical processes. 

Assuming it does, some conclusions may be reached. We begin with only one of the primary processes possessing a 
relevant non-linear term. Its importance for the cyclical growth will depend on its relevance with respect to the coarse 
grained linear approximation of the cyclic propagator. We may conclude that if the linear dominant term originates 
from the same primary process as the non-linearity, the latter will remain relevant. If, however, the non-linear process 
yields the non-dominant part of the free propagator, the non- linearity may be rendered irrelevant, in which case the 
cyclical behavior will be that of the other (linear) process. 

If both primary processes have a non-linear term, a simple behavior is reached if one of them is rendered irrele- 
vant. Then the other one will bequeath its scaling exponents to the cyclical descendant. Theoretically, non- linear 
contributions from both processes may be relevant with various potential outcomes fl2|l . 

To examine these general conclusions, we performed numerical simulations on specific lattice atomistic models in 
2d (preliminary simulations in 3d show a similar behavior). The system size in the simulations was changed between 
128 to 4096 lattice spacings. A typical cycle consisted of a deposition of 5-20 layers (average number of particle 
deposited per site) and desorption of between 10% to 100% of the deposited amount. The maximum number of cycles 
n varied between 500 — 10000, averaged over 50-5000 independent runs, depending on the pairs of primary processes 
and the system size. The growth exponent j3 was extracted for different system size L. The value quoted is from the 
largest L (once (3 became size-independent). From W(L, 00), the saturation width dependence on L, the roughness 
exponent a was derived. In most cases we checked independently the value of a from the scale dependence of the rms 
height-difference. 

For linear primary processes, we looked at the pairwise combinations of RD/EW, its reverse EW/RD, and DT/EW, 
using the standard absorption/desorption algorithms for RD JTJ], EW Q, and DT || Q (see also gj). They all 
showed asymptotic cyclical scaling with EW exponents. This confirms the above conclusions for primary linear pro- 
cesses since EW is the dominating one when paired with RD or DT. We also run the respective adsorption/adsorption 
cycles using the same pairs. The roughness behavior was identical, while their growth velocities were naturally 
different. 

To simulate non-linear processes we included two lattice realizations which belong to the KPZ universality class : 
Ballistic Deposition (BD) (15) |§ Jl6| and the restricted SOS algorithm of Kim and Kosterlitz (KK) @- We have 
generalized both of these algorithms to desorption as well jll[ . Both realizations yielded equivalent behaviors when 
combined with other processes (as they do in simple non-cyclical adsorption or desorption) and we quote the ones 
obtained using the BD (Ballistic Deposition or Desorption) alogarithm. If we combine the non-linear BD (or KK) 
alogarithm with the EW (or DT) linear process to form a cyclical process, the above considerations lead us to expect 
KPZ exponents. Indeed, the KPZ free propagtor is equivalent to the EW (and dominating over the DT) one. The 
exponents obtained are f3 = 0.311(5), a = 0.51(1) for EW/BD , and = 0.322(5), a = 0.50(1) for BD/EW. These 
asymptotic exponents are consistent with the KPZ (3 = 1/3 (and, of course, with a = 1/2, which is the common value 
of EW and KPZ). However, while for BD/EW these values were obtained for all sizes, for EW/BD the exponent f3 
increased slowly with the system size and the effective [3 reached its asymptotic value only for the largest system size 
(L = 4096). Growth/growth cycles yielded results very similar to the growth/recession cycles. 

To look at primary processes with different values of a, a DT («i = 1.5) deposition with ballistic desorption 
(«2 = 0.5) were performed. In Fig. 1 we show the logarithmic dependence of W (roughness) on ln(n), for different 
system sizes L, for this DT/BD process. The inset depicts the logarithmic dependence of Wm (the maximal value of 
the roughness) vs InL. From the graphs we obtain the asymptotic values of the exponents for DT/KPZ: (3 = 0.311(15) 
and a — 0.48(2), both consistent with the KPZ values. 

Finally, simulations of BD/BD and KK/KK (note that they are not time-reversed images of each other because 
of the non-linearity) gave surfaces with KPZ scaling for T\ 7^ T2. For T\ = T2, however, EW behavior was found. 
This follows from the non-linear KPZ terms in the primary processes having the same magnitude but opposite signs. 



Hence they exactly cancel each other in the coarse grained growth equation yielding a EW one. 

Experiments of cyclical growth were performed by metal electrodeposition/dissolution of silver. Electrodeposition 
has been used in recent years to study the scaling behavior of surface growth (see pl| [19J [g(| an d references 
therein) . To explore cyclical growth, multiple electrodeposition/dissolution cycles were carried out on initially vapor- 
deposited silver substrates, ranging from 1 to 20 cycles. The plating solution contained 0.092 M AgBr (silver bromide), 
0.23M(N 114)28203 (ammonium thiosulfate), and 0.17 M (NH4)2SOz (ammonium sulfite). Each cycle consisted of 
plating for 5 minutes followed by 2.5 minutes of clcctrodissolution with a current density of 0.8 mA/cm 2 . Up to 20 
cycles have been performed and the roughness was examined by AFM after n full cycles and after the deposition part 
of the cycles (namely after n+p cycles with p = 2/3) [Ell. Fig. 2 shows a log — log plot of the rms roughness vs. the 
cycle number. The data is consistent with a power-law scaling and the fit yields (3 = 0.48(5). For comparison (3 = 0.71 
for electrodeposition only under the same conditions pl| . Future experimental measurements (on the cyclical as well 
as on the primary processes) will allow more detailed scaling analysis and quantitative comparison of the scaling 
relations with theoretical predictions. 

To summarize, the results of complementary studies of cyclical growth processes were presented, and show them to 
be amenable to scaling analysis. The scaling description holds, provided the time variable is replaced by the number 
of cycles. This conclusion is supported by the initial experimental findings. We have derived the cyclical behavior of 
two alternating linear processes and outlined how the RG approach may be applied in presence of non-linear effects in 
the primary processes. In all systems we have studied analytically or numerically, the related exponents are unaffected 
by the the cycle period T or the relative durations (p and 1-p) of the two processes. One crucial question is how this 
behavior might be affected if the duration of the deposition (and/or desorption) phases are not uniform. We plan to 
address irregular intermittent growth/recession processes in the future |ll| . 

One of us (SR) is grateful to A. Kundagrami for his assistance in the simulations. This research was supported by 
the NSF CMS-9872103 (JJ, SR, YS) and by the Eastman Kodak Company (DGF). 
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FIG. 1. InW (roughness) of the simulated cyclical DT/BD process vs Inn (number of cycles) for different system sizes L 
(inset: InWm (maximal roughness) vs lnL). 

FIG. 2. The roughness vs number of cycles n in the electrochemical cyclical growth of Silver (log-log plot). 
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